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ANALYSIS OF DYNAMIC CAPACITY OF LOW- CONTACT- RATIO SPUR GEARS 

USING LUNDBERG-PALMGREN THEORY 

by John J. Coy, Dennis P. Townsend, and Erwin V. Zaretsky 

Lewis Research Center and 
U. S. Army Air Mobility R&D Laboratory 


SUMMARY 


A concise mathematical model for surface fatigue life of low-contact-ratio spur 
gears was developed. From the model the expected life of the pinion, gear, or gear set 
may be calculated in terms of pinion revolutions. In addition, an equation for the dy- 
namic capacity of the gear set was derived. The dynamic capacity is defined as the 
transmitted tangential load carried by the gear set which will cause the 10-percent fa- 


tigue life to be one million pinion revolutions. 

Experimental gear life data are described for spur gears operating at a Hertz stress 
of 1. 71x10^ newtons per square meter (248 000 psi) and 10 000 revolutions per minute. 
The material was vacuum arc remelted (VAR) AISI 9310 heat treated to a Rockwell C 
hardness of 62. 


The theoretical life was compared to the experimentally determined life. Good 
agreement was obtained when using material constants obtained from life tests of rolling 
element bearings. A previously determined Weibull slope for gears was used to predict 
the survival distribution of the current set of gear specimens. 

The dynamic capacity of the spur gear mesh is given by 


w tM = B 1 cos <n N i 



^jsp (c+h " 1)/2 fte-h - 1 )/ 2 


-1/w 


and the life L corresponding to a particular transmitted load W t is given by 



where is the material constant, cp the pressure angle, N 1 the number of teeth on 
pinion, N 2 the number of teeth on gear, e the Weibull exponent, l j the involute profile 
arc -length on pinion for the zone of single -tooth contact, Sp the curvature sum for the 
point of highest contact stress, f the face width in contact, and p, c, h, and w are the 
material constants. 


INTRODUCTION 


Gears used in aircraft and other applications may fail from scoring, tooth fracture 
due to bending fatigue, or surface pitting fatigue. The scoring type failure is usually 
lubrication related and can be corrected by proper lubricant selection and/or changes in 
gear operating conditions (ref. 1). Tooth breakage is caused by tooth loads which pro- 
duce bending stresses above the endurance limit of the material (ref. 2). It is usually 
accepted that the endurance limit, if it does exist, can be predicted from available 
stress-life (S-N) curves for the material being used (ref. 3). Current methods (ref. 4) 
of predicting gear surface pitting failures are similar to those used for predicting the 
bending fatigue limit. According to the method of reference 4, the maximum surface 
contact stress (Hertz stress) should be limited to a value less than the surface endurance 
limit of the gear material. It is commonly believed that the gears would then have an 
infinite surface pitting life. But based on gearing (refs. 5 and 6) and rolling-element 
bearing life studies (ref. 7), there is no real evidence to support the concept of a surface 
fatigue limit under normal operating conditions of bearings and gears. Rather, it ap- 
pears that all gears even if designed properly to avoid failure by scoring and tooth bend- 
ing fatigue will eventually succumb to surface pitting in much the same way as rolling- 
element bearings. Even so, as Seireg (ref. 8) points out, the maximum Hertzian stress 
is a significant parameter that should be kept within bounds if the gear mesh is to have a 
long enough life. Huffaker (ref. 5) reported that gear life as determined from dynamom- 
eter testing of automotive transmissions is dependent on contact stress according to the 
proportionality 

Life ~ ( i ) X 

\Hertz stress/ 

The failure criterion used was "excessive noise" generated due to pitted gear teeth. 

For stress levels ranging from 1. 4xl0 9 to 2. 6xl0 9 newtons per square meter (200 000 to 
380 000 lb/in. ) the stress-life exponent was found to be 6. 75. The failure distribution 
with stress level was reported to be very similar to that of roller bearings with no evi- 
dence of an endurance limit (i. e. , a stress level for which the gears would have an in- 
finite life). Results of gear testing by Shilke (ref. 6) show similar results for stress 
levels from 1. 4x10 to 1. 9xlO y newtons per square meter (210 000 to 280 000 lb/in. 2 ) 
surface contact stress. Two failure criteria were used: loss of drive and first visible 
pitting under a X10 magnification. Both methods gave similar results. In addition to 
finding that "the pitting fatigue mode of failure appears to have no endurance limit, " 

Shilke concluded that the stress -life exponent is a function of reliability level. 

Reference 9 presents a method of determining the allowable stress in gear teeth for 
infinite life that is consistent with a given level of reliability. However, the method 
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depends on the concept of a surface fatigue endurance limit and is not consistent with the 
findings in references 5 and 6. 

A recent departure in predicting gear surface (pitting) fatigue life from those of 
other investigators was that of Rumbarger (ref. 10). Recognizing the similarity in fail- 
ure between rolling-element bearings and gears, Rumbarger developed a failure model 
(ref. 10) based on the Lundberg-Palmgren life model (ref. 11) for rolling- element bear- 
ings. The Rumbarger model (ref. 10), while a good approach, may have some serious 
limitations as a design tool. In order to apply the model to life predictions several nu- 
merical evaluations of integrals must be carried out. In addition, there is some ques- 
tion regarding the accuracy of the equation as it pertains to the gear tooth profile incor- 
porated in the model, and no full scale gear tests were run to verify the accuracy of the 
model. 

In view of the aforementioned, it becomes the objective of the research reported 
herein (1) to provide a simplified theory for gear surface (pitting) fatigue failure from 
which calculations may readily be made to provide life estimates of low contact ratio 
spur gears and (2) to compare the analytical life prediction with experimental gear life 
data. The method of analysis is based on the rolling- element fatigue theory contained in 
reference 11. Simplifications are incorporated into the failure theory for gears based 
on the observations reported in reference 12, which states that fatigue spalls on gears 
occurred in the region of highest Hertz stress. Generally, no fatigue failures were ob- 
served outside the region of single tooth contact. The results of the theory were com- 
pared to experimental gear fatigue results obtained with gears manufactured from AISI 
9310 material run at 10 000 rpm, a maximum Hertz stress of 1. 71x10^ newtons per 
square meter (248 000 psi), and a temperature of 350 K (170° F). 


LUNDBERG-PALMGREN THEORY 


The fatigue-life model proposed in 1947 by Lundberg (ref. 11) is the commonly ac- 
cepted theory to determine the fatigue life of rolling-element bearings. The probability 
of survival as a function of time or stress cycles is expressed as 




( 1 ) 


where 


S probability of survival 

V volume representation of stress concentration 
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77 millions of stress cycles 

e Weibull exponent 

h, c material dependent exponents 
z Q depth to critical stress 

t q critical stress 

(All symbols are defined in appendix A. ) Hence, if the probability of survival is speci- 
fied, the life 77 for the required reliability can be considered a function of the stressed 
volume V, the maximum critical stress t , and the depth to the critical shearing 
stress z Q . As a result, the proportionality can be written as 


V 


h/e 


T c/e y l/e 


( 2 ) 


In applying this relation to gears in contact, the subscript 1 denotes the pinion and 
2 denotes the gear. The necessary relations pertaining to gear geometry are given in 
appendix B, so that the continuity of the derivation is not interrupted. 


FATIGUE LIFE OF SINGLE PINION TOOTH 

As an expression for the stressed volume in a low contact ratio pinion (1 < M < 2), 
we chose 


v ~ az 0 * 1 (3) 

where li is the involute length of the pinion across the heavy load zone. The choice of 
this length is a simplification of Rumbarger’s (ref. 10) use of the entire involute length 
that passes through the pressure zone. This is a justifiable simplification since opera- 
tion under lighter loads greatly diminishes the probability of failure. Another reason 
for choosing the heavily loaded zone is that all observed failures on the NASA gear tests 
occurred in that zone (mostly where the Hertz load on the pinion is greatest) (ref. 12). 

It is also assumed that the most severe Hertz stress generated in the contact of the teeth 
is acting over the entire load zone of single tooth contact. 

Equation (1) then becomes 
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log — 
S 


c e 7 

T & al i 


( 4 ) 


Using equations (C7), (C9), and (CIO) from Hertz theory (see appendix C) gives 

log i ~ ( 5) 

S \2;rab/ ^ b ^h-l 

Since each tooth is stressed only once per revolution, tj may represent millions of rev- 
olutions of the pinion L. According to appendix B the pinion involute arc length in the 
heavy load zone is given by 


l 1 = A l 


1 2 v 


U1 ' y Ll. 


(6) 


In order that the constants of proportionality and exponents that were found by Lundberg 
and Palmgren (ref. 11) may be used, the following derivations will be closely compared 
to their work. 

Rearranging equation (5) and dropping the constant (3/2 tj) c give 


l T c Q c L e ^ 1 Axc-l/jvc+h-l 
S r h-l \a/ \b/ 


In order to use the same proportionality constant later introduced by Lundberg and called 
the "material constant, " the should be divided by n; the relation is now written as 





( 8 ) 


This equation is valid for both line and point contact. Using the special relations of ap- 
pendix C ((C5), (C7), and (C12)) for line contact gives the 
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This equation agrees in form with equation (46) of reference 11 if note is taken that 


D (c+h-l)/2 D (c-h-l)/2 
a a I 



. (c-h+l)/2 



( 10 ) 


After assuming the probability of survival S = Sj = 0. 9 and designating L = L< as the 
corresponding life for a 90-percent survival probability, the equation becomes 1 


Q (c-h+l)/2 L e 


*1 / 4 \(c-l-h)/ 2 / ff E o 2 p \( c+h - 1 )/2 T c 
n \3f / 1 6 / 


1-1 


h-1 


(ID 


or 


> since T q and are constant for line contact, 


Q (c-h+l)/2 L e _j^l s ^(c+h-l)/2j f ( c -l-h)/2 


(12) 


Next, after defining a material constant Bj and rearranging exponents the form is 
QL 2e/(c-h + l) = Bj p 2p (c+h-l)/2j V(C h+1> f ( c -h-l)/(c-h+l) 

This equation is consistent with equation (48) of reference 11 where the appropriate defi- 
nitions are made for f and D n ; also, because of equation (10) the D is dropped 
Finally, a 


QL^/P) = 


■ t \ -2/(c-h+l) 


~ ( c +h “ 1)/ (c -h+1) f (c-h-l)/(c-h+l) 


(14) 


where 


P = 


c - h + 1 
2e 


(15) 


The life of a single pinion tooth with a 90-percent probability of survival may be de- 
termined by using equation (13) and rearranging terms as follows: 
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L 


1 



S p-(c+h-l)/2e f (c-h-l)/2e 



(16) 


Specific values of the material constant and exponents are not substituted into the 
derivation here in order to maintain generality in the results. 


FATIGUE LIFE OF SINGLE GEAR TOOTH 


If the appropriate variables are substituted into the life equation for the pinion tooth 
(eq. (16)), an equation valid for the gear tooth is obtained: 


L 


2 



Sp"(c + h-l)/2e f (c-h-l)/2e 


Q- p 


(17) 


where the life of the gear tooth is written in terms of gear revolutions, and is the 
involute length of the gear during the single tooth load region. 

The life of the gear tooth in terms of pinion rotations L 2 p may be related to gear 
tooth life L 2 by 



(18) 


Also for the gear tooth as was done for the pinion tooth in equation (6) 


l 


2 





By geometry it can easily be shown that 



(19) 


(20) 


Therefore, 
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l 2 = ~ l l 

N 2 


From equations (16), (17), (18), and (22) the life equation for a single gear tooth 
terms of pinion rotations is 


*“’54 


NAd + e)/e 


where it is assumed that the same load Q is acting on both members. 


LIFE FOR PINION ALONE 

According to the relation for survival probability of a single pinion tooth under coi 
stant service conditions 

log - ~ L e / 

S v 

For a 90-percent survival rate S = Sj = 0. 90 and L = Lj. If the constant of propor- 
tionality K is used, then 


log — = KL, 

HQ A 


Therefore, for any general survival rate we write 


S \L 


io g i = /i\ e log _L 


From basic probability theory, the probability of survival for Nj teeth is the product of 

survival probabilities for the individual teeth. Then the probability of survival for the 
pinion is 


8 


(27) 


S 


P 


N 


= S 


1 


and 


log -L = log ( 1 \ 1 = N, log - = N< M 6 log j- 
S p W 1 S 0. £ 


(28) 


For a 90-percent survival for the pinion, Sp = 0. 9, L = Lp, and 

i / l pY 1 

log — = Njl — — I log-i- 
0.9 \ L l/ 0.9 


(29) 


Then the life for the pinion with a 90-percent survival probability is related to the life of 
a single tooth by 



(30) 


LIFE FOR GEAR ALONE 

Again from the product law, the survival probability for the gear is 



(31) 


where the relation for a single tooth is 

logi= (—] log— (32) 

S \ l 2 ) 0.9 

Then for the gear 

log — = N 2 log/"-\ = N 2 /^l\ log — (33) 

S G W \ L 2/ °- 9 

For a 90-percent survival of the gear Sq = 0. 9 and L = Lq where Lq is the number 
of gear rotations endured: 
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Note that if the life of the gear and gear tooth is given in terms of pinion rotations, then 



where L Gp is the gear life in terms of pinion rotations and L 2p is the life of one gear 
tooth in terms of pinion rotations. 


LIFE OF GEAR SET IN MESH 

From probability theory the probability of survival of the two gears in mesh is 
given by 


S M " S p ' S G 


(36) 


and therefore 


J M, 


J GP/ 


(37) 


For a 90-percent probability of survival for the mesh after L M million revolutions of 
the pinion, substituting equations (30) and (35) into equation (37) gives 


= N, 


+ N„ 


(38) 




TP/ 


J 2P/ 


In this equation it should be noticed that all lives are expressed in the same time base, 
that is, in terms of pinion rotations. 


DYNAMIC CAPACITY OF PINION TOOTH 

The dynamic capacity of a single tooth in the pinion is defined as the transmitted 
tangential load W^ p that may be carried for one million pinion revolutions with a 


10 


90-percent probability of survival. From equation (14), 


Q c i = B il j 


L \-2/(c-h+l) 


2p -(c+h-l)/(c-h+l) f (c-h-l)/(c-h+l) 


W tP = Q cl cos ^ 


(39) 

(40) 


where W^p is the dynamic capacity of a single pinion tooth. 


DYNAMIC CAPACITY OF GEAR TOOTH 

The dynamic capacity of a single gear tooth is defined as the transmitted tangential 
load that may be carried for one million pinion revolutions with a 90 -percent 

probability of survival for that gear tooth. From equations (14) and (23) the gear tooth 
life is 


l -n 

2 p 'Vh, 


(l+e)/e 


B Ai 


l \ - 2/(c-h+l) 


Sp -(c + h-l)/(c-h + l) f (c-h-l)/(c-h+l) q-1 


P 


(41) 


Setting Lgp = 1 and solving equation (41) for the dynamic capacity give 
/N \( 1+e )/ e P fi \-2/(c-h+l) 


QC2 " Ul/ 




Sp -(c+h-l)/(c-h+l) f (c-h-l)/(c-h+l) ^ 


and 


W tG = %2 cos * 

If equations (15), (39), and (40) are used, the dynamic capacity of a single gear 
tooth may be related to the dynamic capacity of a pinion tooth by 


(43) 


W 


tG 


/M 

\ N J 


2(l+e)/(c-h+l) 


W 


tP 


(44) 
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DYNAMIC CAPACITY OF MESH 


In the following derivation the dynamic capacity of the mesh W tM is defined as the 
transmitted tangential load that may be carried for one million pinion revolutions with 
a 90-percent probability of survival of both pinion and gear. The probability of the mesh 
surviving is given by 


N 1 N ? 

S M “ S P ' S G “ S 1 S 2 


(45) 


From equation (9) for a given set of gears the survival probability of the single pinion 
tooth is given by 


where 


log -L = k.Q w L e 
S 1 


(46) 


w = 


c - h + 1 
2 


For the single gear tooth, 


log — = k„Q w L e 
So ^ 


The probability that the mesh will survive is found from equation (45): 


log J_ = (N 1 k 1 + N 2 k 2 )Q w L e 
S M 

In equation (46) if L = 1 and Sj = 0. 9 then Q = Q^. As a result, 


log — = kiOgi 
0.9 1 ^ 1 


In equation (47) if L = 1 and S 2 = 0. 9 then Q = Q c2 _ Thus, 


log — 
0. 9 



(47) 


(48) 
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In equation (48) if L = 1 and S M = 0. 9 then Q - <^, M - Thus 


log — = - k 2 N 2 )Qc M 

0.9 

From the previous three relations 



( 49 ) 


The dynamic capacity of the mesh is then 

w tM - «CM cos V 


(50) 


From equations (15). (39), (42), (46), (49), and (50), the dynamic capacity of the mesh 
or the transmitted tangential load that may be carried for 10 pinion revolutions wi 
90-percent reliability is given by 


f - 


W tM ~ B i cos n N 1 


i+(- 


^LP (C+h " 1)/2 f _(c " h ' 1)/2 


■1/w 


(51) 


A proportion Is set up using equations (48) and (50) whereby the life corresponding 
to any transmitted tangential load may be calculated. The relation is 

iP 


f w. 


L = 


tM’ 


(52) 


W. 


t 


MATERIAL CONSTANTS AND EXPONENTS 

From reference 11 the material constant B was determined from test data on 
Swedish air-melt bearing steel (61. 7 to 64. 5 Rockwell C hardness). This material is 

similar to AISI 52100. The material constant B was determined to be 60 when using 
the gravitational metric system of units (kg and mm) and a 90-percent probability of 
survival. From appendix D, the material constant B x used in this report was e er- 

mined to be 117 when using the gravitational metric system of units (kg an ^mm an 
^percent probability of survival. The B, is 4. 08xl0 8 when using the Internationa 
System (SI) of Units and 102 000 when using the U. S. Customary Units in equation (16). 
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The exponents e, c, h, and p were also determined from test data in reference 11 
or roller bearings with line contact as e = 3/2, c = lo|, h = 2* and p = 3. In a later 
paper by Lundberg and Palmgren (ref. 13), the material constant for line contact was 

as 7 9 ° 8 c = 7° nentS f ° r Hne C ° ntaCt " r ° Uer bearin ^ s were reported 

as e ^ y/a, c = 10£, h = 2|, and p = 4. 

A useful method for adjusting the materiaf constant to account for increases in life 

Thu/ Ts S ‘ eelS ai ’ d be “ er Iubricatlon techniques is presented in reference 14 

atot’fromTt'n 8 Z 3 ‘“l adiUStment factor instead <* changing the material con- 
stant from that originally used in reference 13 would change equations (16) and (39) as 


(l \ _1 /e 

Lj = j Sp -(c+h-l)/2e f (c-h+l)/2e 


Q" 


1 / /^ 1 \" 2 /( c_h+1 ) 

0^1 = D /P BjM 2p -(c+h-l)/( c -h+l) f (c-h-l)/( c -h+l) 


(53) 


(54) 


APPARATUS, SPECIMENS, AND PROCEDURE 

Gear Test Apparatus 

The gear fatigue tests were performed in the NASA Lewis Research Center's gear 
test apparatus (fig. 1). This test rig uses the four-square principle of applying the test 
gear load so that the input drive need only overcome the frictional losses in the system 

„ . ” , C * he teS ‘ r ‘ g ‘ S Sh ° wn ln figure 2 ' 011 Pressure and leakage flow are 

supplied to the load vanes through a shaft seal. As the oil pressure is increased on the 

load vanes inside the slave gear, torque is applied to the shaft. This torque is trans- 
mi e rough the test gears back to the slave gear where an equal but opposite torque 
.s mamiained by the oil pressure. This torque on the test gears, which depends onThe 
hydraulic pressure applied to the load vanes, loads the gear teeth to the desired stress 

/Id l, l I° “ ‘ eSt g6arS Ca " be S ' arted U " der no load ; the load can be applied 
gradually without changing the running track on the gear teeth. 

Separate lubrication systems are provided for the test gears and the main gearbox 

The two lubricant systems are separated at the gearbox shafts by pressurized labyrinth 

sea s. Nitrogen was the seal gas. The test gear lubricant is filtered through a 5-micron 

nomina fiberglass filter. The test lubricant can be heated electrically with an immer- 

Sion eater. The skin temperature of the heater is controlled to prevent overheating the 
test lubricant. 
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A vibration transducer mounted on the gearbox is used to automatically shut off the 
test rig when a gear-surface fatigue occurs. The gearbox is also automatically shut off 
if there is a loss of oil flow to either the main gear box or the test gears, if the test 
gear oil overheats, or if there is a loss of seal gas pressurization. 

The test rig is belt driven and can be operated at several fixed speeds by changing 
pulleys. The operating speed for the tests reported herein was 10 000 rpm. 


Test Lubricant 


All tests were conducted with a single batch of superrefined naphthenic mineral oil 
lubricant having proprietary additives (antiwear, antioxidant, and antifoam). The phys- 
ical properties of this lubricant are summarized in table I. Five percent of an extreme 
pressure additive, designated Anglamol 81 (partial chemical analysis given in table H), 
was added to the lubricant. The lubricant flow rate was held constant at 800 cubic centi- 
meters per minute, and the lubrication was supplied to the inlet mesh of the gear set by 
jet lubrication. The lubricant inlet temperature was constant at 319±6 K (115 ±10 F), 

and the lubricant outlet temperature was nearly constant at 350±3 K (170 ±5 F). This 


outlet temperature was measured at the outlet of the test-gear cover. A nitrogen cover 
gas was used throughout the test as a baseline condition which allowed testing at the 
same conditions at much higher temperatures without oil degradation. This cover gas 
also reduced the effect of the oil additives on the gear surface boundary lubrication by 
reducing the chemical reactivity of the additive-metal system by excluding oxygen 
(ref. 15). 


Test Gears 

Test gears were manufactured from vacuum arc remelted (VAR) AISI 9310 case 
carburized steel to an effective case depth of 1 millimeter (0. 040 in. ). The material 
chemical composition is given in table m and the heat treatment schedule is given in ta- 
ble IV. The nominal Rockwell C hardnesses of the case and core were 62 and 45, re- 
spectively. This material is a commonly used steel in gear manufacture. 

Photomicrographs of the microstructure of the AISI 9310 are shown in figure 3. 
Figure 3(a) shows the high-carbon fine-grained martensitic structure of the hardened 
case of the gear. Figure 3(b) shows the core region of the gear with its softer low- 
carbon refined austenitic grain structure. 

Dimensions for the test gears are given in table V. All gears have a nominal^sur- 
face finish on the tooth face of 0. 406 micrometer (16 pin. ) rms and a standard 20 in- 
volute tooth profile. 


15 


Test Procedure 


The test gears were cleaned to remove the preservative and then assembled on the 
test rig. The test gears were run in an offset condition with a 0. 30-centimeter 
(0. 120-in. ) tooth-surface overlap to give a load surface on the gear face of 0. 28 centi- 
meter (0. 110 m. ) of the 0. 635-centimeter- (0. 250-in. -) wide gear, thereby allowing for 
an edge radius of the gear teeth. By testing both faces of the gears, a total of four fa- 
tigue tests could be run for each set of gears. All tests were run at a load of 1157 new- 
tons per centimeter (661 lb/in. ) for 1 hour. The load was then increased to 5784 new- 
tons per centimeter (3305 Ib/in. ) with a 1. 71xl0 9 newton per square meter (248 000 psi) 
pi c line Hertz stress. At the pitch-line load the tooth bending stress was 24 8xl0 8 
newtons per square meter (35 100 psi) if plain bending is assumed. However, because 
there is an offset load there is an additional stress imposed on the tooth bending stress 
Combining the bending and torsional moments gives a maximum stress of 26 7xl0 8 new- 
tons per square meter (38 700 psi). This bending stress does not consider the effects of 
ip relief which will also increase the bending stress. 

The test gears were operated at 10 000 rpm, which gave a pitch-line velocity of 
n 55 meters P er second (9163 ft/min). A lubricant was supplied to the inlet mesh at 
0 cubic centimeters per minute (0. 21 gal/min) at 319±6 K (115°±10° F). The tests 
were continued 24 hours a day until they were shut down automatically by the vibration- 
detection transducer located on the gearbox, adjacent to the test gears. The lubricant 

?^nn CU l ated thr ° Ugh a 5_miCr ° n fiber S lass filter to remove wear particles. A total 
of 3800 cubic centimeters (1 gal) of lubricant was used for each test; the lubricant was 

discarded, along with the filter element, after each test. Inlet and outlet oil tempera- 
tures were recorded continuously on a strip-chart recorder. 

The pitch-line elastohydrodynamic (EHD) film thickness was calculated by the 
method of reference 16. It was assumed, for this film thickness calculation, that the 
gear temperature at the pitch line was equal to the outlet oil temperature and that the 
mlet oil temperature to the contact zone was equal to the gear temperature, even though 
the oil inlet temperature was considerably lower. It is probable that the gear surface 
temperature could be even higher than the oil outlet temperature, especially at the end 
points of sliding contact. The EHD film thickness for these conditions was computed to 

be 0. 65 micrometer (26 pin. ), which gave a ratio of film thickness to composite sur- 
face roughness (h/cr) of 1. 13. 


RESULTS AND DISCUSSION 

Gear fatigue tests were conducted with gears made from vacuum arc remelt (VAR) 
AISI 9310 steel. Test conditions were a load of 5784 newtons per centimeter (3305 
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lb/in. ), which produced a maximum Hertz stress at the pitch line of 1. 71x10 newtons 
per square meter (248 000 psi), a test speed of 10 000 rpm, and a gear temperature of 
350 K (170° F). A superrefined naphthenic mineral oil was the lubricant. Failure of 
the gears occurred due to surface fatigue pitting. Test results were statistically eval- 
uated using the methods of reference 17. The results of these tests are plotted on 
Weibull coordinates in figure 4. Weibull coordinates are the log-log of the reciproca 
of the probability of survival graduated as the statistical percent of specimens failed 
(ordinate) against the log of time to failure or system life (abscissa). The experimental 
10-percent life or the lift at a 90-percent probability of survival was 11.4 million revo- 

lutions or 19 hours of operation. 

The theoretical 10-percent life for this set of conditions was calculated using equa- 
tion (42) and the exponents h = 2^, c = 10^, e = 3, and constant Bj - 102 000. The 
calculated 10-percent life was 32 million revolutions or 53 hours (see appendix E) where 
h, c, and B x are based on rolling element bearing experience and the Weibull slope e 

is based on gear tests reported in references 5, 6, and 12. 

It should be stated here that the Weibull slope e was assumed to be independent of 

the stress level and reliability level S in the original work referenced (ref. 11). There 

is some evidence suggested in reference 6 showing that the exponent e is dependent on 
the stress level. However, the value of e used previously is representative of data by 
Shilke (ref. 6) and other NASA tests where the stress level was the same as in the tests 

reported here. 

The predicted life can be considered a reasonably good engineering approximation 
to the experimental life results. However, the theoretical prediction does not consider 
material and processing factors such as material type, melting practice, or heat 
treatment - nor does it consider environmental factors such as lubrication and tempera- 
ture. All these factors are known to be extremely important in their effect on rolling- 
element bearing life (ref. 14). There is no reason why these effects to determine gear 
life should be significantly different from those used to determine bearing life. More 
test data obtained with gear specimens under various test conditions and different ma- 
terials and lubricants are required to establish and/or affirm the material constant Bj 
and the exponents c, h, and e for gears. However, the results presented herein sup- 
port the use of the statistical methods presented for predicting spur gear fatigue life 
with a standard involute profile. 


SUMMARY OF RESULTS 

An analytical model was developed to determine the fatigue life and dynamic capac- 
ity of low contact ratio spur gears. The analytical results were compared with expen- 
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mental gear life data obtained with a group of vacuum arc remelted (VAR) AISI 9310 
spur gears. The test gears had a standard 20° involute profile and a 8. 89 -centimeter 
(3. 5-in. ) pitch diameter. Test conditions were a maximum Hertz stress of 1. 71X10 9 
newtons per square meter (248 000 psi), a speed of 10 000 rpm, and a temperature of 
350 K (170 F). The lubricant was a superrefined naphthenic mineral oil with an addi- 
tive package. The following results were obtained: 

1. There was good agreement between the predicted gear mesh life and the experi- 
mental life results. - F 

2. The experimentally determined Weibull slope e for a sample of spur gears and 
t e material constant Bj and exponents h and c from roller bearing life tests were 
used successfully to predict gear life. However, further experimental work is needed 

o give statistical significance to those exponents and the material constant. 

3. The dynamic capacity of the spur gear mesh is given by 


W, 


tM “ ®1 cos 


■ /n \ e *i/z \ y i/w 

l+( — - j |i S/) (c+h-l)/2 f -(c-h-l)/2 ( 

WJw J 


and the life corresponding to a particular transmitted load is given by 



where 



B i 

material constant 

= 

<p 

pressure angle 


N i 

number of teeth on pinion 

=- 

N 2 

number of teeth on gear 

r 

e 

Weibull exponent 


h 

involute profile arc-length on pinion for zone of single-tooth contact 


Sp 

curvature sum for point of highest contact stress 
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f 


face width in contact 


p, c, h, w material constants 

Lewis Research Center, 

National Aeronautics and Space Administration, 
and 

U. S. Army Air Mobility R&D Laboratory, 
Cleveland, Ohio, March 26, 1975, 

505-04. 


APPENDIX A 


SYMBOLS 

half of major axis of Hertzian contact, m (in. ) 

material constant defined by eq. (116) of ref. 11 

material constant defined by eq. (13) 

half of minor axis of Hertzian contact, m (in. ) 

gear center distance, m (in. ) 

orthogonal shear stress exponent 

life improvement factor defined by eq. (53) 

rolling element diameter, m (in. ) 

raceway diameter (in direction of rolling), m (in. ) 

pitch diameter of roll bodies 

Young's modulus, N/m (psi) 

defined by eq. (C3) 

Weibull exponent 
face width of tooth, m (in. ) 
depth to critical stress exponent 
proportionality constants 
life in millions of revolutions 
involute profile arc length, m (in. ) 
contact ratio 
number of teeth 

diametral pitch, teeth/m (teeth/in. ) 
load-life exponent defined by eq. (15) 
base pitch, m/tooth (in. /tooth) 
circumferential pitch, m/tooth (in. /tooth) 
load normal to involute profile, N (lb) 
dynamic capacity for normal loading, N (lb) 
pitch circle radius, m (in. ) 


r„ addendum circle radius, m (in.) 

a. 

r^ base circle radius, m (in. ) 

S probability of survival 

T defined in eq. (Cll) 

U number of contact cycles per bearing revolution 

V volume, m^ (in. j 

transmitted tangential load, N (lb) 
w defined by eq. (50) 

Z contact path length, m (in.) 

z Q depth of occurrence of critical shearing stress, m (in.) 

o' angle of approach, rad 

a’ angle of recess, rad 

a contact angle, rad 

/3 contact roll angle, rad 

y total contact roll angle, rad 

5 precontact roll angle, rad 

Z number of roll bodies per row 

£ defined in eq. (CIO) 

rj millions of stress cycles 

9j^ base circle roll angle at entry to single tooth load zone, rad 

base circle roll angle at exit from single tooth load zone, rad 
pinion base circle roll angle, rad 
p, v defined in eq. (C2) 
p curvature radius, m (in. ) 

Dp curvature sum, m"^ (in. "*) 

9 

t q maximum subsurface orthogonal reversing shear stress, N/m (psi) 

cp pressure angle, rad 

Subscripts: 

G gear ' 

H high load 
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L low load 

M mesh of pinion and gear 

P pinion 

1 reference to driving member, pinion 

2 reference to driven gear 
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APPENDIX B 


GEAR GEOMETRY 

CALCULATION OF PINION INVOLUTE ARC LENGTH AS 
FUNCTION OF PINION ROLL ANGLE 

The radius of curvature at point A is r^j©^ the arc length along the involute is 

/ 9 2^/2 

l , and the differential element of arc length d£ is (dx + dy J . The coordinates of a 
typical point on the involute are (x, y) for the roll angle 8^ (fig. 5): 


: = r bl 

sin 

V 

- tjA cos 0j 

= r bl 

sin 

V 

- r^^ cos 0 1 

H 

X 

II 

(sin 

V 

- 0J cos 0j) 

r bl + 

r bl 

cos 

01 + tjA sin 0j 

r bl + 

r bl 

cos 

0j + r^i^j sin 0j 


(Bl) 


= r^^cos + 8 1 sin 0j - 1) 
dx = r^j^cos d0j - d0j cos 0j + 0^ sin 8 j d0j) 


(B2) 


= sin d0j 


dy = r b ^(-sin 8 j d0^ + d0^ sin 0j + 0^ cos 0^ d0j) 


= cos 0j d0j 


dZ = V^ r bl 0 l sin 0 1 + ( r bl 0 l cos 0 1 “ r bl 0 l 


(B3) 


After integrating between any two angular positions of the pinion, the increment of in- 
volute length is 
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Al 


1 “ 





(B4) 


Calculation of Base Pitch 


The base pitch is defined as the distance between teeth along the pressure line. 
Base pitch is also the distance from a point on one tooth to the corresponding point on 
the next tooth measured along the base circle (see fig. 6): 


Base pitch = p^ = 


2 ^ r bi 



2^ r bj cos <p 


N 


1 



= P c cos (P 


(B5) 


Contact Ratio 

The path of contact Z is the distance between where the two addendum circles 
cross the pressure line as shown in figure 7: 

Addendum radius = r + — 

P 


Similarly, 


Z = eg 


1 \ 2 
etr, = 'll I r 2 + — j - (r 2 cos <p) 


ft, = r 2 sin <p 


1 2 

ef = (r 2 + — 1 - (r 2 cos <p) - r 2 sin cp 


1 \<i 2 

fg = + ~ ) - (rj cos <p) - rj sin <p 
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Let 


C = Center distance = r i + r 2 


Z = +ij - <rj cos V ) 2 + y( r 2+ l) 2 - (r 2 cos »>) 2 - C sin ^ (B6) 

The contact ratio M is defined as the average number of gear teeth in contact at one 
time. It is calculated by 



When two teeth are in contact the load is assumed to be equally shared by the teeth. 
For low contact ratio gears (1 < M < 2) there are three load zones. The first occurs as 
the teeth are coming into contact. At that instant there is also another pair of teeth in 
contact on the line of action. When the contact occurs near the pitch point there is only 
one pair of teeth in contact and the total load is carried by a single pair of teeth. The 
third zone is similar to the first and occurs when the teeth come out of contact. The 
zones are shown in figure 8. 

Two sets of teeth are in contact during roll angles /3 L and only one set during ro 
angle In the figures 9 and 10 two teeth on each gear are shown by their contacting 

face profile only. There are two positions shown: with teeth coming into contact in fig- 
ure 9 and teeth going out of contact in figure 10. 

The roll angle for which two teeth are in contact is calculated by 


0L1 



(B7) 


And the heavily loaded zone 0 Hl may be computed by subtracting the two lightly loaded 
zone angles from the total contact angle: 


Z - 2(Z - 

0H1 = r l ' 2 ^L1 = 

r bl 


2p b ~ Z 
r bl 


(B8) 
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Calculation of Angles of Approach and Recess 


The angle of approach a is the angle between a line joining the gear centers and a 
second line which is drawn from the gear center to the intersection of the involute pro- 
. e with the pitch circle. This intersection is chosen for the instant when the gears are 
coming into contact as shown in figure 11. Angles of approach may be defined for each 
gear; in general, a 1 ±a 2 - 

As the pinion rotates through angle a x pbint e will advance to point f; that is 
the contact point will move along the line of action: 


ef « et 2 - ft 2 . 


1 \ 2 
,r 2 + ~ ) " ( r 2 cos <P) ~ r 2 sin <p 


Also since the line of action can be thought of as a string being wound onto the base 
circles as the gears turn in mesh. 


ef = r 


bl^l 


The previous two equations may be solved to 
gear: 


give the approach angles for the pinion and 


a i = 



r 2 + ~ ) ‘ ( r 2 cos ^) 2 - r 2 sin cp 


l bl 


ato = 


' l a l 


(B9) 


The angle of recess a ' is defined as the angle included between the line of centers 

t e line from the gear center to the intersection of the involute profile with the pitch 
circle at the instant when the gears are going out of mesh as shown in figure 12. 

The rotation angle of the pinion will wrap the length fi of the line of action 
onto the base circle so 


r bi 
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'min' miii I in I 11111 ||| mi, 


But 


fg = 


*i f = 




2 

- (r^ cos cp) 


2 


- sin cp 


Therefore, 


and 




cos <p) 2 - rj sin <p 
r bl 


'N 


> 



The total angle of contact is 


yj = a?! + a\ 

Finally, 6 X is defined as the angle of roll that will bring the pinion tooth from 
point to a contact with the mating gear tooth (fig. 13): 


6 


1 



But 


et l = *1*2 ” e *2 


where 


*1*2 = r l Sin + r 2 sin ^ = C s * n ^ 


(BIO) 


the base 
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and 


et 2 = 


= [f 2 + p) " (r2 C ° S 


1/2 


Therefore, 



2 + “ ) ~ ( r 2 cos <PY 


(Bll) 


Also note that 


h f 

- 6 j + a j = tan <p (B12) 

r bl 


which is evident from figure 13. The required angles defining the length of the most 

heavily loaded length of the involute profile may be written by referring to figure 8. The 
needed equations are 


0 L1 “ 6 1 + / 3 L1 

(B13) 

0 U1 = d Ll + ^H1 

(B14) 
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APPENDIX C 


STRESS EQUATIONS FOR GEARS IN LINE CONTACT 
The curvature sum in terms of pinion roll angle (fig. 14) is determined first: 



_ 1 _ 

p 2 


where 


Sp = 


Pi = r hl e 1 


P 2 = *1*2 ~Pi = Csmcp- r^j 

1 C sin <p - r bl e l + r bl g l . 


C sin <p 


r bl 0 l C sin <p - r fel 0j r bl 0 1 (C sin cp - r^j) r bl 0j(C sin (p - 


(Cl) 


It was determined in reference 11 that 



where a and b are the major and minor half axes of the contact ellipse and 


E o = 


m 2 E 


m 


(C2) 


(C3) 


where 1/m is Poisson’s ratio. 

Thomas and Hoersch (ref. 18) have shown that when line contact is approached in the 
limit (a - 00 , Q ~ °°) the general equations for point contact may be extended for that 
case. With appropriate changes in notation from their equation (84) it is found that 
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(C4) 


2n~b 2 a rt> 2 f 
3Q 2 Q 

From this it follows that for a uniform distribution of pressure on the face of the gear 
the semimajor axis may be taken as 


a = — f (C5) 

4 

Also, from reference 18 it can be shown that 

lim p, 2 = 2 (C6) 

b/a-0 jr 

By using equations (C2), (C5), and (C6) we calculate b as 

b = (C7) 

\7rfE 0 Sp 

Equations (C5) and (C7) are used in the development of equation (9) in the report. 

For the general case of point contact the maximum contact stress occurs at the 
ellipse center and is given by 


q = —0- (C8) 

2irab 

where Q is the normal load applied to the body. For the line contact situation in gears 
the maximum Hertz stress on the surface is 

q=^ (C9) 

nfb 

According to reference 11 the cause of fatigue flaking of bearings is the maximum re- 
versing orthogonal shear stress that occurs at a depth z Q below the surface and has an 
amplitude that varies between ±t q where 


IS 

o 

II 

tr 

(CIO) 

r 0 =Tq 

(Cll) 
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In general, $ and T depend on the shape of the contact ellipse and for the special case 
of line contact 



a 


T Q = 0. 2500 > 

£ 0 = 0. 5000 


(C12) 


If nonferrous gear materials are used, then the contact dimensions a and b and the 
stress field parameters T and £ will change due to the change in E Q for the different 
material. 
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APPENDIX D 


MATERIAL CONSTANT 

The tests run by Lundberg and Palmgren for SKF bearing steel (AISI 52100) with 
Rockwell C hardnesses of 61. 7 to 64. 5 yield a material constant of 


B = 60 (kg-mm units) 


for the dynamic capacity of line contact given by equation (116) of reference 11: 

2/9 

d m ) \ d m/ 


BIT 


D cos *\ 29 / 27 /D 

X d 29/27 z 7/9 z -1/3 




/ D cos a 
f l ±_S 


™\l/3 


(Dl) 




Equation (Dl) is valid for a roller bearing. The task now is to find the value of B ^ to 
be used in the dynamic capacity equation for gears which will be applicable for U. S. 
Customary (lb and in. ) units. In order to simplify the referencing of equations from the 
work of Lundberg and Palmgren all numbered equations in this appendix refer to the 
numerations used in reference 11. According to equation (58) of reference 11, 

/n \ 2 / 9 

Q c = B l( D aS p)- 35 /27^\ u -1/3 d 29/27,7/9 (D2 

Then with equation (107) 

u = — fl ± — cos 
2 \ 

and equation (106) where F(b/a) = 1 for line contact (b/a = 0) gives 



° a 2p 


2 

D cos a 
1 t 


(D4) 
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Using equations (D3) and (D4) in equation (D2) gives 


Qc 


= B 



-35/27 


D„ cos a 


m 


'D 


D_ 


2/9 


g / ± D a cos a 

2 \ d m ' 


-1/3 


d 29/27,7/9 (w) 


a a 


From equations (103) and (104) 


D = d 1 ± 
n ml d 


D„ cos of' 
a 


m 


(D6) 


After substituting equation (D6) into (D5) the result is 


Q„ = 
^c 


[B, 


,1/3 


'It 


n 29/27 

D„ cos a\ 


m 


^D„ V 


2/9 


j3 5/27 j i Q cos q,\ 1/3 \^m/ 


D 29/ 2 7 i 7/9 g -V3 (D7) 

1 a a 


1 ± 


m 


Therefore, 


B 


1 


2 1 / 3 

2 35/27 


60 


B, = 116. 95 

This concludes the determination of B^. Equation (D2) may be written as 

2/9 

Q^(kg) = 116. 95(D a Sp)~ 35 / 27 ^^ u" 1 / 3 [D a (mm)] 29 / 27 [Z a (mm)] 7 / 9 (kg) (D9) 

If conversion factors are used in equation (D9), the value of Bj which is valid for 
U. S. Customary Units (lbf and in. ) may be determined: 
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n 2/9 

O (lb) = — — x 116. 95 X (D Sp)' 35 / 27 u _1 /3 

0. 4536 kg a D 


n 


x jiyin.) X - 2 M mm j // p a (in. ) x 2 - 5 -^ mm j ?/9 (kg) ( D i 0 ) 


Therefore, 


B 1 = :S 50/27 - 102 000 (U. S. Customary (lb-in. ) Units) (Dll) 

Similarly, Bj = 4. 08xl0 8 when SI units (N-m) are used. In a later publication (ref. 13) 
Lundberg and Palmgren report a material constant B of 56. 2. When this is reduced as 
the value of 60 was before, then the material constant Bj = 95 500 (U. S. Customary 
(lb-in.) Units) or 3. 82xl0 8 (SI units (N-m)) results. 
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APPENDIX E 


CALCULATION OF BASIC DYNAMIC CAPACITY AND LIFE OF TEST GEARS 


Gear data: 


Nj = N 2 = 28 

f = 0.00279 m (0. 11 in.) 

C = 0. 08890 m (3. 500 in. ) 

P = 8 

r l = r 2 = 0. 04445 m (1. 750 in. ) 


Test load: 


W t = 1615 N (363 lb) 


Material constants: 


Bj = 102 000 

c = 10— 

3 

h = 2^ 


e = 3 


From figure 7, 


■ 2 = r bl = rj cos <p = 0.04177 m (1.6444 in.) 


By equation (Bll), 


§1 = 0. 180186 rad 



By equation (B5), 


Pb = 0- 009373 m (0. 369016 in. ) 


By equation (B6), 


Z = 0.015353 m (0.604450 in.) 


By equation (B7), 


By equation (B13), 


By equation (B8), 


*L1 = °- 143168 rad 


Or i = 0. 323354 rad 


= 0. 081232 rad 


By equation (B14), 


- 0. 404585 rad 

The worst Hertzian stress is considered to occur at the lowest point of single tooth con 

tact on the pinion. Therefore, the curvature sum is calculated for that point, and the 

stress is assumed to be constant while only one pair of teeth are in contact. By equa- 
tion (Cl), 


2p = 133.2 m' 1 (3.383627 in." 1 ) 


By equation (6), 


By equation (46), 


l j = 0. 001235 m (0. 048620 in. ) 


w = 4. 5 
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By equation (51), 


W tM - 16284 N (3660. 81 lb) 


By equation (15), 


p = 1. 5 


By equation (52), 


L = 32. 03 million revolutions of the pinion 
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TABLE I. - PROPERTIES OF SUPERREFINED, NAPHTHENIC, 
MINERAL -OIL TEST LUBRICANT 


2 

Kinematic viscosity, cm /sec (cS) at 


266 K (20° F) 


311 K (100° F) . . 

7^y 1 n’2 

372 K (210° F) .... 

7 7x1 n ^ ( 7 7\ 

477 K (400° F) . . 

i r v i n “ 2 / 1 a \ 

Flash point, K (°F) . . AOO om 

Autoignition temperature, K (°F) 

♦ . . . . -r vc/ 

664 (735) 

Pour point, K (°F). 


Density at 289 K (60° F), g/cm 3 

0 RftQQ 

Vapor pressure at 311 K (100° F), mm Hg (or torr) 

....... 0.01 

Thermal conductivity at 311 K (100° F), J/(m)(sec)(K) (Btu/(hr)(ft)(°F)) . 

. . 0. 04 (0. 0725) 

Specific heat at 311 K (100° F), J/(kg)(K) (Btu/(lb)(°F)) 

... 582 (0. 450) 


TABLE II. - PROPERTIES OF LUBRICANT ADDITIVE ANGLAMOL 81 


Percent phosphorous by weight 0. 66 

Percent sulfur by weight 13. 41 

Specific gravity 0. 982 

Kinematic viscosity at 372 K (210° F), cm 2 /sec (cS) ... 29. 5xl0' 2 (29. 5) 


TABLE m. - CHEMICAL COMPOSITION OF VAR AISI 9310 
GEAR MATERIALS BY PERCENT WEIGHT 


Element 

C 

Mn 



Si 

Ni 

Cr 

Mo 

Cu 

P and S 

Weight 

percent 

0. 10 

0. 63 

0.27 

3. 22 

1. 21 

0. 12 

0. 13 

0. 005 
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TABLE IV. - HEAT TREATMENT PROCESS FOR VACUUM 


ARC REMELTED (VAR) AISI 9310 


Step 

Process 

Temperature 

Time, 

hr 

K 

°F 

1 

Carburize 

1172 

1650 

8 

2 

Air cool to room temperature 

— 

— 

— 

3 

Copper plate all over 

— 

— 

— 

4 

Reheat 

922 

1200 

2. 5 

5 

Air cool to room temperature 

— 

— 

— 

6 

Austenitize 

1117 

1550 

2. 5 

7 

Oil quench 

— 

— 

— 

8 

Subzero cool 

189 

-120 

3. 5 

9 

Double temperature 

450 

350 

2 each 

10 

Finish grind 

— 

— 

— 

11 

Stress relieve 

450 

350 

2 


TABLE V. - SPUR GEAR DATA 
[Gear tolerance per ASMA class 12.] 


. 28 

. 8 


Number of teeth 

Diametral pitch . * 

Circular pitch, cm (in.) 0.9975 (0.3927) 

Whole depth, cm (in. ) 0. 762 (0. 300) 

Addendum, cm (in.) 0.318 (0. 125) 

Chordal tooth thickness reference, cm (in.) 0. 485 (0. 191) 

20 

Pressure angle, deg 

Pitch diameter, cm (in. ) 8 - 890 50 °) 

Outside diameter, cm (in. ) 9. 525 (3. 750) 

Root fillet, cm (in. ) 0. 102 to 0. 152 (0. 04 to 0. 06) 

Measurement over pins, cm (in. ) ... 9. 603 to 9. 630 (3. 7807 to 3. 7915) 

Pin diameter, cm (in. ) 549 

Backlash reference, cm (in. ) 0.0254 (0.010) 

Tip relief, cm (in. ) 0. 001 to 0. 0015 (0.0004 to 0. 0006) 















Figure 7. - Length of contact path. 









